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Abstract 


An  instrument  has  been  developed  to  measure  the  turbulent  variations 
of  velocity  and  temperature  in  a river,  and  with  it  some  observations  have 
been  made  in  the  Kennebec  Estuary  of  the  vertical  transport  of  heat  and 
momentum.  A comparison  is  made  with  the  values  determined  from  the  aver- 
age distribution  of  these  properties  by  the  method  of  Jacobsen.  Only  a 


small  part  of  the  turbulence  transports  momentum  and  heat;  and  the  rela- 


tion between  the  ratio 


/^h  / ^ 


v-  and  the  Richardson  number  is 


found  to  be  similar  to  that  obtained  by  Jacobsen  and  Taylor. 


Introduction 


The  mixing  of  fresh  and  salt  water  has  previously  been  measured  in 
estuaries  by  Jacobsen  (1918)  and  Pritchard  (1952)  from  the  salinity  and 
velocity  distribution  averaged  over  several  tidal  cycles.  Pritchard's 
method  is  more  refined,  but  Jacobsen's  method  illustrates  the  principle 
well.  It  is  assumed  that,  when  averaged  over  a tidal  cycle,  the  only  im- 
portant terns  in  the  horizontal  equation  of  motion  are  the  horizontal 
pressure  gradient  and  the  shearing  stresses,  and  that  the  vertical  equa- 
tion is  simply  the  hydrostatic  equation.  Then  cross-differentiation 
eliminates  the  pressure  term  and 

d JT  J_  If 
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where  ^ is  the  horizontal  stress  along  the  estuary,  2 is  measured 
vertically  upward  from  the  bottom,  Q is  the  acceleration  due  to  gravity 
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and  is  density.  The  assumptions  amount  to  using  steady  flow  equations 
for  conditions  that  are  actually  periodic  and  can  only  be  averaged  over 
an  integral  number  of  periods.  The  stress  is  obtained  by  direct  numeri- 
cal integration  with  respect  to  It . Assumptions  also  have  to  be  made 
about  boundary  conditions  which  enter  as  constants  of  integration.  Simi- 
larly, the  vertical  flux  of  salt  (or  heat)  is  obtained  by  integration 


T 
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where  ^ is  salinity  (or  heat  content).  If  observations  of  the  time  mean 
vertical  gradients  of  velocity  and  salinity  are  also  available  the  eddy 
coefficient  of  viscosity  and  $ of  diffusivity  may  be  computed: 

^ = z /(dWd*)  ; ^ --  Y/(*l>*) 

Jacobsen  (1918)  carried  out  such  computations,  and  Taylor  (1931)  showed 
that  the  ratio  x is  related  to  the  Richardson  number 
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Because  these  confutations  are  a rather  indirect  way  of  measuring 
the  vertical  turbulent  transports  of  properties,  it  is  appropriate  to  make 
measurements  of  the  turbulent  fluctuations  themselves  at  various  places 
and  times  in  an  estuary  to  compute  the  turbulent  transfers  directly  from 
the  formulae  ~C  =■  — j?  K.  ' {Ar  ' 


r=  f T~  er 


and 
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where  the  primes  denote  deviations  from  the  mean,  U/^  is  vertical  velocity, 
T is  temperature,  and  Cp,  is  specific  heat.  A strictly  complete  com- 
parison of  these  direct  measurements  with  those  obtained  by  the  Jacobsen 
method  is  of  course  impossible  because  a single  ship  can  observe  only  at 
one  place  at  a time,  whereas  the  Jacobsen  method  averages  over  both  space 
and  time.  The  limited  comparison  possible,  however,  is  interesting  because 
it  permits  an  inspection,  for  the  first  time,  of  the  size,  intensity,  and 
period  of  the  very  eddy  processes  doing  the  transport. 

General  Description  of  the  Kennebec  Estuary 

Three  short  cruises  to  the  Kennebec  Estuary  were  made  at  the  sugges- 
tion of  Dr.  A.  C.  Redfield  who  suggested  that  on  account  of  the  geometri- 
cal simplicity  of  its  shape,  it  might  be  an  ideal  location  for  various 
types  of  experiments.  The  valley  is  straight  and  of  fairly  uniform  width 
and  depth  (Figure  l);  only  two  sharp  bends  occur  (Doubling  Bends);  and 
there  is  one  very  marked  transition  in  width  (at  Bluff  Head).  Autumn  dis- 

•3  T 

charge  is  about  290  nr  sec  ; a tide  of  3 m at  the  mouth  produces  currents 
of  up  to  2 m sec“^  in  the  estuary.  These  produce  considerable  mixing. 
During  the  first  short  cruise,  which  was  of  an  exploratory  nature,  some 
idea  of  the  velocities  and  densities  of  the  estuary  were  obtained.  Two 
stations  were  occupied  for  a tidal  cycle:  one  at  the  mouth,  the  other 

at  Bath,  about  20  km  upstream.  Velocity  was  measured  with  a Roberts 
meter,  density  by  hydrometer.  Figure  2 shows  the  main  results:  slight 

stratification,  currents  strong  enough  for  easy  propeller  measurement. 

The  estuary  is  somewhat  more  mixed  vertically  than  Randers fjord,  but 
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clearly  should  offer  a good  place  to  make  vertical  turbulent  transport 
measurements. 


Second  Cruise 

A second  cruise  was  undertaken  (R/V  ASTERIA8)  to  test  an  instrument 
devised  for  measuring  turbulent  fluctuations  in  horizontal  velocity  and 
in  inclination  of  flow.  A propeller  meter  (von  Arx,  1950)  of  the  induc- 
tion-coil type  was  used.  The  pulsed  voltage  output  was  rectified  and 
smoothed  by  a condenser  and  then  recorded  on  a Speedomax  recorder.  In- 
clination of  the  flow  to  the  vertical  was  measured  by  a vane  rotating  on 
a horizontal  axis  ahead  of  its  leading  edge.  The  vane  was  attached  to  a 
potentiometer,  the  output  of  which  was  recorded  on  another  Speedomax. 

The  direction  of  flow  in  a horizontal  plane  was  not  measured,  and  because 
the  ocean  happened  to  be  at  the  same  temperature  as  the  fresh  water,  no 
temperature  fluctuations  were  present. 

Various  types  of  suspension  were  tried:  (i)  suspension  by  its  flexi- 

ble cable,  (ii)  fixed  to  a tripod  resting  on  the  bottom,  (iii)  fixed  to  a 
vertical  mast  projecting  down  into  the  water  from  the  ship.  The  records 
obtained  by  methods  (i)  and  (iii)  show  ary  rolling  or  pitching  of  the 
ship.  The  advantage  of  knowing  the  exact  depth  and  position  of  the  in- 
strument in  method  (iii)  was  offset  by  the  difficult  problem  of  rigging 
it  in  a strong  current.  Method  (ii)  is  confined  to  measurements  at  the 
bottom.  Therefore  method  (i)  was  generally  employed.  In  calm  weather 
there  does  not  appear  to  be  any  perceptible  difference  in  records  ob- 
tained by  the  three  different  methods. 
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After  many  instrumental  difficulties,  one  velocity  profile  was  obtained 

at  Station  5C  (depth  15  m)  in  perfectly  calm  weather  on  an  ebb  tide.  A ten- 

minute  record  was  taken  at  20  cm  depth,  3 m,  6 m,  and  9 m,  the  meter  being 

suspended  on  its  cable.  At  depths  below  9 m there  was  so  little  velocity 

that  the  meter  did  not  register;  a salinity  reading  showed  that  this  was  an 

example  of  a salt  layer  underlying  the  fast  ebbing  fresher  water. 

The  records  of  horizontal  velocity  to  and  inclination  of  flow  were  • 

analyzed  as  follows.  A two-minute  period  was  selected  and  the  starting 

and  end  points  identified  in  both  records.  Ordinates  were  scaled  off  each 

at  three-second  intervals  and  the  mean  values  U~  and  ^ computed.  The 

variations  ( v /=  - tc)  and  {•&'-  &■  - &)  were  taken  for  each  ordi- 

/ / 

nate,  and  the  product  ^ M.  O'  formed,  which  is  proportional  to  the 
, / 

M (***  ( ur  = instantaneous  vertical  velocity)  if  it  is  assumed  that  the 

mean  flow  is  horizontal  (e.g.  &■  = 90°  from  vertical)  and  that  -Q-  is  small 

/ / / o 

so  that  sin  & = O . In  fact,  G never  exceeded  10  , and  there  is  no 

evidence  to  suppose  that  the  mean  flow  was  inclined  to  the  horizontal. 

The  mean  value  of  n 'u'  over  the  U0  ordinates  was  converted  to 

the  Reynold’s  stress  = — j*  "v*  using  the  suitable  constants 

for  the  velocity  and  angle  scales.  The  computations  were  repeated  twice 

for  succeeding  two-minute  intervals  and  the  values  of  7"  are  tabulated  in 

Table  I. 

The  differences  of  "£  from  the  successive  periods  of  computation 
show  presumably  how  some  of  the  momentum  is  transported  vertically  by  ed- 
dies of  scale  too  large  to  be  observed  in  the  two-minute  period.  However, 
the  order  of  the  stresses  is  about  what  might  be  expected  as  the  wall 
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First 

period 


Second 

period 


Third 

period 


Table  I 

Station  5C  - Second  Cruise 
Values  of  Shear  Stress 


Depth  cm 


20 


4A. 

cm  sec-'1' 

12li 

T 

cfynes  cm-^ 

6.7 

cm  sec-"1 

127 

V 

dynes  cm“^ 

22.1 

lZ 

cm  sec-'1 

120 

nr 

dynes  cm-<^ 

lU.l 

cm  sec-1 

12U 

nr 

dynes  cm- 

111. 2 

120 

550 

730 

98 

62 

33 

16.1 

8.3 

3.2 

107 

51 

37 

21.5 

32.9 

12.2 

100 

67 

la 

9.3 

8.2 

12.6 

102 

60 

37 

15.6 

16.U 

9.3 

Mean 
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stress  of  a stream  800  cm  deep  flowing  over  a very  smooth  surface  at  120  cm 
sec”^.  A point  of  some  interest  lies  in  the  high  stresses  observed  quite 
near  the  surface.  It  is  clear  that  the  stress  at  the  very  surface  must  be 
zero  on  this  windless  day;  yet  at  only  20  cm  below  it  was  as  high  as  at  all 
other  points  through  the  fluid.  Either  the  instrument  must  in  some  subtle 
way  give  a false  reading,  possibly  due  to  its  large  size  compared  with  its 
submergence,  or  some  hitherto  unsuspected  complication  is  caused  by  the 
free  surface. 

During  this  cruise,  a survey  was  made  to  find  the  horizontal  salinity 
gradients  along  the  estuary.  It  was  carried  out  on  an  ebb  tide,  starting 
at  the  furthest  upstream  station  at  U2  hours  after  high  water  at  the  mouth, 
and  arriving  at  the  mouth  of  the  estuary  § hour  before  low  water  slack  tide 
and  6|  hours  after  high  water.  The  return  was  made  on  the  flood  tide.  The 
data  are  presented  in  Figure  3,  the  conductivity  of  the  water  (mho  cm“^) 
being  plotted  against  distance  downstream,  for  the  water  at  both  the  sur- 
face and  9 m below.  The  sudden  increase  of  conductivity  at  Station  2 is 
due  to  salt  water  being  discharged  into  the  main  river  from  the  tributary 
on  the  left  bank,  which  rapidly  sinks  to  the  bottom  of  the  river;  other- 
wise the  curves  for  both  depths  are  somewhat  parallel,  and  the  salinity  in- 
creases rather  uniformly  toward  the  sea.  It  is  interesting  to  observe 
that  the  two  sharp  bends  (Doubling  Bends)  and  narrow  portion  downstream 
from  them  (between  Stations  3 and  5)  do  not  appear  to  change  the  salinity 
gradient  very  much.  This  is  surprising,  because  in  this  region  the  appar- 
ent turbulence  is  great,  great  "boils"  of  water  rising  to  the  surface 
giving  characteristic  circular  patches  of  still  water.  Such  turbulence  is 
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absent  elsewhere.  It  might  be  thought  that  these  "boils'*  are  doing  a great 
deal  of  mixing  of  the  heavier  salt  water  below  with  the  fresher  water 
above;  if  in  fact  they  were  doing  so,  the  vertical  salinity  (or  conductivity) 
gradients  would  be  reduced,  and  the  9 m conductivity  would  be  more  nearly 
the  surface  one.  No  such  effect  is  noticeable  and  a possibility  exists 
that  visible  turbulence  in  the  form  of  "boils"  does  not  necessarily  cause 
salt  mixing. 


Third  Cruise  (R/V  CARYN) 

The  object  of  this  cruise  was  to  continue  the  observations  of  turbu- 
lent transfer  processes  of  salt  and  momentum,  and  to  investigate  if  the 
sudden  widening  of  the  river  at  Bluff  Head  acts  as  a control  section  on 
the  velocity  and  density  stratification  upstream  of  it. 

The  velocity  and  vertical  direction  meter  used  in  the  second  cruise 
was  used  again,  with  an  additional  element  to  measure  temperature  fluctua- 
tions, which  was  a "thermistor"  bead  connected  in  a bridge  circuit  to 
another  Speedomax  recorder.  As  the  ocean  water  was  this  time  some  10°F 
warmer  than  the  fresh  river  water,  the  temperature  of  the  mixed  waters  in 
the  estuary  was  a tracer  of  the  salinity.  A correlation,  Figure  U,  between 
the  water  temperature  and  density  was  obtained  from  samples  taken  both  at 
the  surface  and  below  by  a Nansen  bottle.  Density  was  measured  by  hydro- 
meter. 

The  whole  meter  assembly  was  suspended  on  the  cable,  with  a heavy 
"depressor"  weight  of  about  19  kgm  secured  below  from  a davit  which  pro- 
jected some  1.8  m amidships.  There  is,  therefore,  a possibility  that  the 
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meter  was  affected  by  the  proximity  of  the  ship's  side  when  it  was  close  to 
the  surface.  Also,  in  a swift  current,  and  the  meter  well  immersed,  the 
1.2  cm  diameter  cable  vibrates  at  about  a ^ second  period  with  a surface 
amplitude  of  several  cm.  It  is  not  known  how  the  meter  performance  was 
affected.  Some  uncertainty  was  experienced  in  deciding  the  direction  of 
flow  under  the  surface,  particularly  at  Station  $B  where  on  the  ebb  there 
was  separating  flow  from  Bluff  Head.  A large  drag  or  current  cross  sus- 
pended on  a thin  wire  sometimes  was  used  to  show  this  direction  relative 
to  the  ship. 

Observations  were  made  over  approximately  a whole  tidal  cycle  at  each 
of  four  stations  5,  !?A,  $B,  5>C  as  shown  in  the  ma,f  Figure  5.  Stations 
3>B  and  5C  are  immediately  downstream  of  the  enlargement  of  the  river  cross 
section,  is  in  the  narrow  section,  and  5 is  well  downstream.  The  bot- 
tom is  occasionally  rocky  but  the  soundings  show  that  there  are  no  great 
discontinuities,  and  there  is  no  great  change  of  depth  at  Bluff  Head. 

The  observations  were  done  in  windless  weather,  with  the  air  temperature 

/ 

nearly  the  same  as  that  of  the  water,  except  for  those  at  Station  5 when 
there  was  a fresh  breeze  against  the  current  and  therefore  there  were  waves 
on  the  water  surface.  The  ship's  motion  was  just  perceptible  on  deck,  but 
was  very  obvious  on  the  records  of  the  inclination  vane,  which  were  so 
obscured  as  to  make  them  useless  for  the  turbulence  analysis.  At  each 
station  a sounding  was  made  at  two  hourly  intervals,  when  the  combined 
instrument  was  lowered  from  the  surface  in  3 m increments  of  cable  length. 
At  each  increment  a record  was  taken  for  about  ten  minutes,  and  the  cable 
angle  recorded,  in  order  to  know  the  depth  of  the  instrument.  During  the 


"surface"  measurement  (really  with  meter  center  20;  cm  from  surface) -the 
surface  velocity  was  independently  measured  by  timing  a .floating  object, 
passing  down  the  length  of  the  ship,  and  the  s ur f ace'  density  and  tempera- 
ture observed. 

The  velocity  and  temperature  records  were  fir§t  examined  to  find  how 
the  velocity  and  temperature  gradients  changed  during- the  .'tidal  cycle. 

The  mean  velocity  “K  and  temperature  f”  were  estimated  by  eye  • for • every 
ten-minute  measurement,  and  were  plotted  as  velocity  and  temperature  pro- 
files against  depth  in  Figure  6.  Here  the  horizontal  abscissa  is  the 'time 
of  the  middle  of  the  sounding  relative  to  high  water  at  Bluff  Head,.-  and 
the  vertical  lines  in  each  profile  have  been  placed  at  this  time.  All 
four  stations  are  shown  in  the  figure,  and  the  velocities  and  temperatures 
are  all  plotted  at  the  same  scale.  An  uncertainty  sometimes  arises,  and 
is  shown  by  dotted  lines  when  the  surface  velocity  was  repeated  after  a 
sounding  (e.g.  perhaps  1 hour  after  it  was  initially  measured)  and  was 
found  to  have  changed.  This  was,  of  course,  noticeable  near  high  and  low 
water,  and  is  well  illustrated  by  a curious  phenomenon  noticed  at  Station 
5 on  an  ebb  tide.  The  meter  was  at  20  cm  depth  and  registering  about  60 
cm  per  water  speed,  when  a line  of  foam  was  noticed  approaching  the  ship. 
As  the  line  passed  the  meter,  the  velocity  quite  suddenly  increased  some 
5>C$  and  the  temperature  fell,  showing  that  a wedge  of  quiite  differently 
constituted  water  was  passing  (see  Figure  7).  Subsequent  "fronts"  or  sud- 
den increases  of  velocity  were  also  noticed  and  these  render  doubtful  the 
soundings  at  those  times,  e.g.  Station  $ at  +3\  hours;  .Station  5B  at  -2|- 
hours  and  +2  hours;  Station  5C  at  +1^  hours  from  high  water’.-  •’ 
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The  velocity,  pro  files  of  Figure  6 show  how  flood  begins  first  in  the 
salt  lower  water.,,  but  ebb  begins  first  at  the  surface.  This  is  different 
from  the  phase  lag  due  to  friction  alone  (Lamb,  1932)  where  both  flood  and 
ebb. begin  at  the  bottom. 

Analysis  of  Third  Cruise  Data 

The  records  were  further  analyzed  for  their  turbulence  characteris- 
tics in  the  same  manner  as  described  for  the  records  of  the  second  cruise. 
In  this  case  the  heat  transfer  2T  = f*0**  was  computed  from  the 

temperature  record  | and  ov*  = “U  . A great  deal  of  rather  subjec- 
tive selection  was  required  to  determine  suitable  parts  of  the  record  for 
these  techniques.  Records  which  had  produced  mean  profiles  of  velocity 
or  temperature  of  an  S shape  were  rejected,  for  such  profiles  cannot  be 

reliably  used  to  determine  gradients and  ZJ-  with  any  certainty; 

those  taken  when  there  were  surface  waves  present  (and  therefore  the  ship 
was  moving)  were  rejected  because  the  movement  showed  in  the  inclination 
vane  and  velocity  records,  thus  giving  some  uncertainty;  those  which 
showed  no  variation  in  the  temperature  of  the  water  were  rejected,  for 
without  these  fluctuations  there  will  be  no  heat  transfer;  and  records 
Were  rejected  in  which  the  mean  velocity  had  altered  in  the  course  of 
taking  -the  whole  profile,  such  as  that  shown  in  Figure  7.  When  all  these 
conditions  were  applied,  only  four  profiles  were  regarded  as  wholly  suit- 
able, and  three-minute  portions  of  their  records  were  analyzed  for  the 
depths  shown  in  the  accompanying  Table  II.  The  data  of  *U  , and  ~J~ 

were  then  examined  by  a statistical  method^"  to  find  the  probability  that 

, • - - ■ — — - — ■ — ■ ■ ■ - - - — - - 

1The  te$t  involved  is  that  given  by  Fisher  (19U1),  p.  202,  who  gives  tables 
. of  the  probability  of  a certain  correlation  by  a random  arrangement  as  a 
furiictibn  of  the  number- of  observations. 
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Depth  m 


0.20 

2.13 

5.80 

8.20 


0.20 

3.0 

6.1 


3.0 

6.1 


0.20 

U.85 

8.5 


Table  II 

Computed  Turbulence  Data 


oc  cm  sec"^- 

T °F 

T dynes  cm-^ 

Yf  cals 

Significance  Test* 

cm"*  sec"x 

for  x 

for  f 

Station  $k 

2 hrs.  before  HW 

130 

37 .1*7 

6.82 

0.001*9 

4 1/100 

> i/io 

126 

37.63 

25.3 

0.250 

41/100 

ca  1/20 

113 

38.57 

16.7 

0.179 

41/20 

ca  1/20 

116 

38.91 

17.1 

-0.123 

41/50 

ca  1/20 

Station  5C 

5 hrs.  before  HW 

83.9 

37.00 

1.1*5 

0.0138 

>1/10 

66.2 

37.03 

-0.0326 

-0.0026 

>1/10 

>1/10 

30.1* 

37.18 

-1.7U 

+0.0057 

>1/10 

Station  5C 

5 hrs.  after  HW 

115 

37.96 

32.9 

0.102 

4 1/10 

>1/10 

97.9 

37.93 

20.5 

0.0626 

>1/10 

Station  5C 

3^  hrs. 

after  HW 

ll*0 

39.^9 

1*5.3 

0.106? 

\ 

119 

39.61* 

23.6 

0.211  V 

< 1/100 > 

< 1/100 

85.0 

39.70 

62.2 

0.192  J 

J 

‘^Probability  that  T and  'V'  are  simply  result  of  random  correlation 
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the  computed  values  of  "T"  and  IT  could  have  appeared  from  a completely 

/ / — t 

random  set  of  values  of  U , <*r  and  ( . Where  this  probability  was 
better  than  1/10  the  computed  value  of  X”  or  ^ was  considered  to  be  so 
uncertain  as  not  to  warrant  the  drawing  of  deductions  which  are  presented 
in  Table  III.  This  is  tantamount  to  deciding  that  in  the  discarded  cases 
there  was  so  much  turbulence  present  which  did  not  transport  momentum  or 
heat  that  it  masked  the  turbulence  that  was  in  fact  transporting  these 
quantities. 

The  deductions  are  made  in  groups  of  three  depths,  the  velocity  and 
temperature  gradients  being  found  from  the  extremes,  and  the  transport  of 
momentum  and  heat  from  the  middle  depths.  Other  methods  of  analysis  are 
of  course  possible,  which  may  in  the  last  case  reduce  the  anomaly  of  JA  s 
being  larger  than  a situation  which  is  physically  difficult  to  ex- 

plain. For  example,  in  Table  IIIA,  the  same  deductions  are  made,  but  now 
the  values  of  % and  ^ used  to  find  and  are  the  mean  values 

from  all  the  three  depths.  The  Richardson  number  remains  unchanged. 


Kennebec  Estuary  Turbulent  Transfer  Coefficients  by 
Method  of  Jacobsen 


The  data  in  Figures  3,  k,  and  6 may  be  used  to  obtain  average  values 
of  the  turbulent  exchange  coefficients  by  the  method  of  Jacobsen.  As  a 
working  hypothesis  may  be  considered  as  independent  of  . 

/ ax 

In  order  to  take  account  of  the  integration  constants,  it  may  be  assumed 
that  T" = 0 at  i - 0y  the  bottom  and  the  free  surface,  in  which 

case  t is  a maximum  at  mid-depth  and 
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Table  III 

Deduced  Quantities  from  Turbulence  Analysis 


Depth  m 


Station  5A  2 hrs.  before  HW 


830 

127 

1010 

85 

0.15U  0.191 

0.081*  0.0525 


0.20 

li.85 

8.5 


Station  5C  3§  hrs.  after  HW 


356  830 


2.33 


7.08 


J 


Table  IIIA 


Depth  m 


Station  5A  2 hrs.  before  HW 


0.20  j 

2-13 ; 

535 

72.5 

0.136 

5.80  jy 

1190 

U8.5 

0.01^08 

8.20  J 

Station  5C  3|  hrs.  after  HW 


Stability 


0.191 

0.0525 


- 1 6 - 


T = ^ ^ 

"•*  T 9 

This  parabolic  distribution  of  r is  of  course  an  approximation  to  the 
much  more  complicated  distribution  shown  by  the  turbulence  analysis 
actually  to  exist.  From  Figure  3 estimates  were  made  of 


-Lit 

/ a-x 


*=  8.3  x 10~^  cm"^;  h "=  20 


m 


and  from  Figure  6 


tc 

a'z 


-2  -1 

= 3.3  x 10  sec 


Hence  X = h dynes  cm”^  and  ^ M - r/(j^/d-i) 


2 -1 
= 120  cm  sec 


Similarly  the  mean  temperature  gradient  is  dT /d% 

cm”'*',  and  <LT  dbt  - 1/1200  °F  cm”'*'.  The  turbulent  transport  of  heat  at 

~ *T 

mid-depth  is  approximately  *f  = 


= 3.3  x 10-6  °F 


is  the  mean  dis- 


i<r  ax  "here  $*■ 1 

charge  per  unit  width  of  the  lower  upstream  moving  heavy  water.  The  value 
of  can  be  roughly  estimated  from  the  conservation  of  salt.  Inasmuch 

as  the  average  density  of  the  bottom  layer  is  1.010,  and  that  of  the  surface 
layer  about  1.008,  fa  must  be  about  U times  that  of  the  fresh  water  alone, 

^ o . As  ^ 0 is  the  river  discharge  (290  m^  sec”'*')  divided  by  the  aver- 
age width  of  the  channel  ( *=  580  m)}  thus  £4.  = 2 ra^  sec"'*'  = 2 x 10^  cm**  sec"^ 


and  2^  = °»°7  °F  cals 


cm” ^ sec”-*- 


and 


= T 


/£L 

/ 


7ft  2 -1 

= 7o  cm  sec 
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The  above  values  are  now  combined  to  give 

-£j£/ '(jzt)  =1-9 

as  the  mean  mixing  constants  for  the  estuary. 


Discussion  of  Turbulent  Mixing  Computations 


The 


values  of  the  ratio  obtained  by  our  turbulent  transfer 


measurements  have  teen  plotted  logarithmically  against  the  reciprocal  of 
the  Richardson  number  on  the  graph  Figure  8 which  is  based  on  that  given 
by  Anderson  and  Hunk  (19U8) . The  value  found  by  the  Jacobsen  method  for 
the  Kennebec  is  plotted,  as  well  as  the  original  Taylor- Jacobsen  values. 

It  is  clear  from  the  figure  that  a stable  system  tends  to  inhibit  mixing 
and  so  give  a low  j , a result  well  known  qualitatively.  Con- 
versely, if  the  velocity  gradients  are  great  and  density  gradient  small, 
the  heat  (or  salinity)  and  momentum  are  equally  well  mixed  and  Nh  - 
1.0.  The  mechanism  of  mixing  of  salt  cannot  therefore  be  the  same  as  the 
mixing  of  momentum  (Taylor,  1931).  Although  it  is  by  no  means  original 
with  us,  the  following  physical  explanation  seems  quite  reasonable:  that 

a great  mass  of  heavy,  slow,  liquid  is  torn  from  a lower  layer  by  turbulent 
pressure  fluctuations;  it  then  is  ejected  into  the  faster,  higher  liquid 
where  it  speeds  up  and  so  takes  up  momentum  by  pressure  forces  on  it.  Dur- 
ing this  time  the  outer  surface  is  being  frittered  away  and  so  some  mixing 
of  the  mass  is  caused.  However,  in  very  stable  systems  not  all  the  liquid 
is  thus  mixed,  and  the  reduced  mass  falls  back  again  to  its  parent  liquid 
but  now  at  an  increased  speed,  where  its  momentum  is  absorbed.  Inspection 
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of  the  velocity  and  temperature  records  shows  that  periods  of  eddies  from  1 
to  1$0  seconds  could  sometimes  be  seen,  though  it  is  possible  that  the  longer 
periods  may  not  have  been  adequately  averaged  in  the  180  second  period  of 
the  computations. 


The  comparatively  good  agreement  of  our  values  of 


f'/f*' 


with  those 


obtained  by  Jacobsen  is  perhaps  a vindication  of  the  latter  method,  depen- 
dent as  it  is  on  a not  entirely  justifiable  averaging  process  over  a tidal 
cycle.  It  seems  from  the  turbulence  analysis  that  in  part  of  this  estuary 
at  some  states  of  tide  the  mixing  of  salt  and  fresh  water  is  good  and 
Hr  44  — 1.0;  and  that  at  times  the  stratification  is  pronounced,  when 

(see  figure  8).  It  is  gratifying,  therefore,  to  find  that 
the  apparent  mean  value  of  and  Richardson's  number,  found  by 

Jacobsen's  method,  lies  between  the  sounding  analyzed  which  had  been  classi- 
fied as  well  stratified,  and  the  sounding  which  had  weak  stability  and  was 
therefore  well  mixed.  Although  the  numerical  data  for  the  Jacobsen  method 
are  approximate,  the  value  is  not  changed  much  relative  to  those  found  by 
the  turbulence  analysis  if  other  approximations  are  taken;  and  we  make  the 
.very  tentative  conclusion  that  to  a certain  degree  the  mixing  and  stability 
of  an  estuary  are  rather  linear,  and  that  average  values  seem  to  have  some 
significance  relative  to  spot  values. 

It  will  be  seen  from  Table  II  that  high  stresses  are  again  found  near 
the  surface.  Whether  this  effect  is  a false  one,  due  purely  to  the  instru- 
ment's limitations  or  the  method  of  suspension,  remains  to  be  seen.  At 
any  rate,  it  is  another  reason  for  further  study  of  the  surface  layers  of 
an  open  stream. 
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Some  Suggestions  for  Future  Work 

We  think  that  further  study  on  the  mixing  processes  in  estuaries  can 
develop  along  two  distinct  lines. 

1.  Further,  more  detailed  observations  are  required  on  the  gross  fea- 
tures of  the  salinity  and  velocity  patterns  in  estuaries  with  different 
sorts  of  control  sections.  More  refined  methods  are  required  of  determin- 
ing velocity  and  density  gradients  quickly,  so  that  the  essentially  un- 
steady tidal  flow  of  the  estuary  may  be  sampled  at  all  depths  at  nearly 
the  same  time.  This  requirement  is  probably  to  be  met  by  having  several 
meters  registering  at  the  same  time,  at  different  depths,  though  these 
meters  need  only  give  mean  velocities  and  not  the  turbulent  fluctuations. 

2.  The  small  scale  turbulent  details  of  the  flow  require  study  to 
know  more  exactly  the  mixing  of  the  fluids  under  given  mean  velocity  and 
density  gradients.  This  study  involves  measurements  both  close  to  the 
bottom  and  to  the  free  surface,  and  requires  instruments  with  a quick 
enough  response  to  register  the  turbulent  fluctuations  of  the  current. 

It  is  desirable,  again,  to  have  a number  of  instruments  at  different 
depths  operating  simultaneously,  though  this  will  undoubtedly  add  to  the 
complications  of  the  experiment.  We  think  that  difficulties  due  to  non- 
steady conditions  would  be  eliminated  if  the  experiments  were  done  in  a 
non-tidal,  steady  flow  channel,  in  effect  a large  flume,  if  such  an 
estuary  could  be  found. 
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